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(57) Abstract 



The Conway-nSloane decoding method presents a general premise that for any given highly complex block code, an attempt should 
be made to find a subcode of that given block code that is less complex to decode. The given complex block code may then be maximum 
liklihood decoded with reduced complexity by cycling over the subcode and its cosets. Extending this general premise to Reed-Muller 
codes, and in particular to a second order Reed-Muller (or Kerdock) code, reduced maximum likelihood decoding complexity is obtained 
by cycling over a found subcode (and its cosets) of the second order Reed-Muller (or Kerdock) code. A first order Reed-Muller code 
comprises just such a suitable subcode which, like the parity check for the Golay codes analyzed by Conway and Sloane, has a known less 
complex maximum likelihood decoding solution utilizing a fast Hadamard Transform. 
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METHOD AND APPARATUS FOR DECODING 
SECOND ORDER REED-MULLER CODES 

BACKGROUND OF THE INVENTION 
5 Techruq<a Field of tfre In vention 

The present invention relates to the decoding of 
digital data transmitted over a communications channel 
and, in particular, to a low complexity maximum likelihood 
decoder for second order Reed-Muller codes . 

10 Description of Related Art 

There exist many applications where large volumes of 
digital data must be transmitted and received in a 
substantially error free manner. In telecommunications 
systems, in particular, it is imperative that the 

15 reception of digital data be accomplished as reliably as 
is possible. Reliable communication of digital data is 
difficult, however, because the communications channels 
(including radio frequency, fiber optic, coaxial cable, 
and twisted copper wire) utilized for data transmission 

20 are plagued by error introducing factors. For example, 
such errors may be attributable to transient conditions 
in the channel (like interference, noise or multi-path 
fading) . The influence of such factors results in 
instances where the digital data is not transmitted 

25 properly or cannot be reliably received. 

Considerable attention has been directed toward 
overcoming this problem and reducing the number of errors 
incurred when transmitting data. One option involves 
increasing transmitter power. However, this is typically 

30 not practical due to limitations regarding transmitter 
electronics, regulations on peak power transmission, and 
the added expense involved in increasing power levels. 
A preferable alternative option for combating noise on the 
communications channel is to introduce redundancy in the 

35 transmitted message which is used at the receiver to 
correct introduced errors. Such redundancy is typically 
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implemented through the use of error control coding 
(channel codes) . A preferable alternative option for 
combating fading on the communications channel is to use 
an interleaver to reorder the data prior to transmission 

5 over the channel. As a result, many communications 
systems now utilize a combination of error control 
coding/decoding and interleaving/ de - interleaving 

processes to protect against the effects of interference,, 
noise or multi-path fading on the communications channel. 

10 Because of implementation complexity concerns, the 

error control decoder typically used comprises a soft 
decision decoder (and, in particular, an errors and 
erasures decoder) . Such decoders exploit reliability 
values output from a demodulator in estimating the 

15 transmitted codeword. In the absence of fading, and in 
the presence of Gaussian noise, the optimal soft decision 
decoder is the maximum likelihood decoder. It is also 
typically the best decoder in the presence of fading 
(assuming a good estimate of the fading is available) . 

20 For a general block code, however, maximum likelihood 
decoding can be hopelessly complex to implement. 
Accordingly, a need exists for a less complex maximum 
likelihood decoding scheme for implementation in 
connection with soft decision decoding of block codes. 

25 For the special case of the (24,12) Golay code and 

the (23,12) extended Golay code, a maximum likelihood 
decoder having a very low complexity has been devised by 
Conway and Sloane (see, IEEE Trans. Infor. Theory, vol. 
32, pp. 41-50, 1986). The premise behind the Conway - 

30 Sloane decoding method is that for a given Golay code, an 
attempt is made to find a subcode of that given Golay code 
that is easy to decode. The given Golay code may then be 
decoded by cycling, to achieve a lower overall complexity, 
over the subcode and its cosets. For the (24,12) Golay 

35 code, for example, it is noted that there is a subcode 
thereof which is equivalent to a parity check code . Such 
a parity check code presents a trivial decoding challenge. 
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A more complete explanation of the operation of the 
Conway- Sloane decoding method may be obtained by referring 
to the previously mentioned IEEE article, or to U.S. 
Application for Patent Serial No. 08/768,530, filed 

5 December 18, 1996, by Ali S. Khayrallah, et al., the 
disclosures of which are hereby incorporated by reference. 

The disclosed Conway- Sloane decoding method is 
limited on its face to application in connection with the. 
Golay and extended Golay codes. Furthermore, an extension 

10 of the Conway- Sloane decoding method has been proposed by 
the above -referenced U.S. Application for Patent in 
decoding various shortened (19,7), (18,7) and (18,6) Golay 
codes. In spite of these advances relating to more 
efficient decoding of Golay codes, a need still exists for 

15 less complex maximum likelihood decoding schemes 
specifically addressing other types of codes. In 
particular, there is a need for such a scheme in 
connection with the decoding of second order Reed-Muller 
codes . 

20 

SUMMARY OF THE INVENTION 

To address the foregoing and other problems, the 
decoder of the present invention processes a received 
vector of second order Reed-Muller (or Kerdock) encoded 

25 information by cycling over all of the cosets of a subcode 
of the second order Reed-Muller (or Kerdock) code, and 
performing a Fast Hadamard Transformation on the received 
vector with each cycle to obtain an intermediate codeword 
guess with respect to each examined coset . Once cycling 

30 through the cosets is completed, a final codeword guess 
(comprising the best one of the intermediate codeword 
guesses) of the received vector is obtained. The 
originally encoded information bits are then recovered 
from this final codeword guess. 



WO 99/09655 



PCT/US98/16780 



-4- 



BRIEF DESCRIPTION OF THE DRAWINGS 

A more complete understanding of the method and 
apparatus of the present invention may be acquired by 
5 reference to the following Detailed Description when taken 
in conjunction with the accompanying Drawings wherein: 

FIGURE 1 is a functional block diagram of a block 
coding system of the present invention; and 

FIGURE 2 is a flow diagram illustrating the decoding 
10 process implemented by the channel decoder of FIGURE 1 for 
second order Reed-Muller codes. 

DETAILED DESCRIPTION OF THE DRAWINGS 

Reference is now made to FIGURE 1, wherein there is 

15 shown a block diagram of a block coding system 100 of the 
present invention. The system 100 includes a transmit 
side 102 and a receive side 104. On the transmit side 
102, the system 100 includes an (n,k;d) linear block 
channel encoder 106 wherein a block of "k" information 

20 bits received from an information source encoder 108 is 
encoded to output a codeword of "n" bits in length 
(wherein n>k) . The channel encoder 106 preferably 
implements an error control code. An example of the 
information source encoder 108 is a vocoder or data 

25 compressor. The code words output from the channel 
encoder 106 are then rearranged by an interleaver 110. 
A modulator 112 then maps the rearranged code words into 
waveforms suited for transmission over a communications 
channel 114. An example of such a modulator 112 is any 

30 known modulator having an M-ary signal constellation (like 
quadrature amplitude modulation (QAM) or phase shift 
keying (PSK) ) . The communications channel 114 is 
typically a wireless medium which suffers from many error 
and distortion introducing problems including fading, 

35 interference and noise. 

On the receive side 104, the system 100 includes an 
appropriate demodulator 116 that demodulates the 
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communications channel 114 transmitted communication and 
outputs estimates of the rearranged code words. The 
demodulator 116 further outputs soft information 
comprising a reliability value for each estimated code 
5 word. The reliability values are indicative of the level 
of confidence expressed by the demodulator 116 in its 
estimate of a particular received and demodulated code 
word. Demodulators 116 producing code word estimates and 
reliability values are well known in the art, and thus 
10 will not be further described. Examples of such 
demodulators 116 include: a maximum aposterior probability 
(MAP) demodulator, or a soft output Viterbi algorithm 
(SOVA) demodulator/decoder. The estimated code words and 
corresponding soft information are then reordered (i.e., 
15 de- rearranged) by a de- inter leaver 118. An (n,k;d) linear 
block channel decoder 120 then processes the reordered 
estimated code words and corresponding soft information 
to generate estimates of the information bits for output 
to an information source decoder 122. The channel decoder 
20 120 preferably comprises a maximum likelihood decoder for 
the selected error control code which utilizes soft 
decision decoding. 

In accordance with the present invention, the channel 
encoder 106 implements a second order Reed-Muller code. 
25 An (n,k;d) Reed-Muller code is determined by two 
parameters: m and r^m, wherein: 
n = 2 m 
k = 1 + 
d = 2 m " r 

30 For r=l, a first order Reed-Muller code is presented, some 
examples of which include the (8,4;4) Reed-Muller code, 
the (16,5;8) Reed-Muller code and the (32,6;16) Reed- 
Muller code. For r=2, on the other hand, a second order 
Reed-Muller code is presented, some examples of which 

35 include the (16, 11; 4) Reed-Muller code, and the (32, 16; 8) 
Reed-Muller code. The channel encoder 106 of FIGURE 1 
preferably, but not exclusively, implements as its second 
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order code the (32,16;8) Reed-Muller code. As an 
alternative , the channel encoder 106 may implement a class 
of nonlinear codes (commonly referred to as Kerdock codes) 
sandwiched between first and second order Reed-Muller 
5 codes . 

Although the present invention focuses on the 
decoding of second order Reed-Muller (or Kerdock) codes, 
a more complete understanding of the present invention may 
be obtained from further understanding a well known less 
10 complex procedure used for maximum likelihood decoding of 
first order Reed-Muller codes. A Sylvester-type Hadamard 
matrix may be constructed as follows: 

H 0 =[l] 

H , H , 

m-1 m-1 

H = 

m H -H , 

m-1 m-l_ 



(i) 
(2) 



The columns of H (k) are the code words of the Hadamard code 
in a specific order. One way to construct a first order 

15 (2 ffl / m+l / 2 m ~ 1 ) Reed-Muller code is to take as code words the 

rows of H m (referred to as y 0 , . . . ,y 2 ra -i) , and their 2 m 
additive inverses (referred to as y 2 m , . . . ,Y?™ ¥l -\) . Note 
here that code word y t corresponds to the information bits 
making up the natural binary representation of i. If one 

20 of these code words is generated and transmitted, it is 
received as a receive vector r= (r 0 , . . . , r n . x ) . The maximum 

likelihood decoding of the received vector reduces to the 
following process: 

- compute the fast Hadamard Transform (HT) 
25 r = (r ft , . . . , r J = r H , 

- find I = arg max |r| , 

- if f > 0, then guess y x was the sent code word, 

- otherwise, guess y I+2 m was sent the sent code word, 
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- recover the information bits from I and J of the 
guessed code word. 
Maximum likelihood decoding of a first order Reed-Muller 
code using the above less complex process requires 
5 m x (2 m - 1) additions and 2 ra - 1 comparisons. This is 
much more efficient than the "brute force" (i.e., complex) 
maximum likelihood decoding process which requires 
2 m«i x (2 m - i) additions and 2 m - 1 comparisons. 

Reference is now once again made to the reduced 

10 complexity maximum likelihood decoder for Golay codes 
devised by Conway and Sloane (see, IEEE Trans. Infor. 
Theory, vol. 32, pp. 41-50, 1986). It will be remembered 
that in accordance with the disclosed Conway-Sloane 
procedure an attempt is made to find a subcode of a given 

15 Golay code that is less complex to decode. The given 
Golay code may then be maximum likelihood decoded by 
cycling, to achieve a lower overall complexity, over the 
found subcode and its cosets. The foregoing Conway-Sloane 
decoding method presents a general premise that for any 

20 given highly complex block code, an attempt should be made 
to find a subcode of that given block code that is less 
complex to decode. The given complex block code may then 
be maximum likelihood decoded with reduced complexity by 
cycling over the subcode and its cosets. Extending this 

25 general premise to Reed-Muller codes, and in particular 
to the second order Reed-Muller code used by the channel 
encoder 106, it is proposed that reduced maximum 
likelihood decoding complexity may be obtained by cycling 
over a found subcode (and its cosets) of the second order 

30 Reed-Muller code. The first order Reed-Muller code 
(discussed above) comprises just such a suitable subcode 
which, like the parity check for the Golay codes analyzed 
by Conway and Sloane, and in accordance with the previous 
decoding description, has a known less complex maximum 

35 likelihood decoding solution utilizing a fast Hadamard 
Transform (as discussed above) . 
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Consideration is now specifically given to a second 
order (n,k;d) Reed-Muller code C wherein: 
n = 2 m 

k=:l+m + m(m - 1) 
2 

5 d = 2 m " r 

This second order Reed-Muller code C contains a first 
order (n,k (1, ;d 1 ) Reed-Muller subcode C a) wherein: 
n = 2 m 
k (l) = 1 + m 
10 d 1 = 2 m " r 

The code words of the second order Reed-Muller code C may 
be generated as the union of the code words of the first 
order Reed-Muller subcode C a> and its cosets in C. To 
produce the cosets, consideration must be given to another 
15 (n,k <2) ;d 2 ) Reed-Muller subcode C (2) wherein: 

n = 2 m 

k (2) = k _ k (i) = m(m-l) 

2 2 
d 2 = 2 m ~ r 

The code words of the Reed-Muller subcode C t2) are denoted 
20 z Q , • • • ,z <a> , with z 0 being the all zero code word in 

binary representation. The cosets of the Reed-Muller 
subcode C U1 are then defined as follows: 
- for each z 3 in the subcode C (2) , 

C d,j) = + y ^ e C (l)j 

25 These cosets C u ' j> are disjoint, and their union is equal 
to the second order Reed-Muller code C. 

With the subcode comprising the first order Reed- 
Muller subcode C U) and its cosets C a ' j \ the general 
premise underlying the Conway-Sloane decoding method may 

30 be extended to decoding the second order Reed-Muller code 
C. If a code word in the second order Reed-Muller code 
C is generated and transmitted, it is received as a 
receive vector r= (r 0/ . . . ,r H ) .. Instead of comparing r 

directly to a coset C (1,j) , a sign change is first performed 
35 to obtain r': 
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r' = r * — j (3) 



wherein: the multiplication (*) is performed on an 
element by element basis; and 
ij comprises in bipolar form (where 0 in 

binary becomes +1, and 1 in binary becomes 
5 -1) . 

It is noted here that a code word yi + Zj in C, with 
component y A in C (1) and z jin C (2) , corresponds to k 
information bits where the first k m bits comprise the 
natural binary representation of i, and the remaining k <2) 
10 bits comprise the natural binary representation of j . The 
maximum likelihood decoding of the received vector then 
reduces to the following process: 

- 1=0, p=0 

- for j=0 : 2 k - 1 

15 - change signs using Equation (3) , 

- compute the fast Hadamard Transform 

- find i' = arg max |r| , 

- if \r I /\ > p , then 
20 - p = |* x ,|, 

- if f z > 0 then 1=1', otherwise I=I'+2 m , 
~ J«j, 

- end 

- end 

25 - guess y x + Zj was the sent code word, 

- recover the information bits from I and J where the 
first k U) bits comprise the natural binary representation 
of I, and the remaining k (2) bits comprise the natural 
binary representation of J. 

30 Maximum likelihood decoding of a second order Reed-Muller 
code using the above less complex process requires 
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2 k<2> x m x 2 m additions and 2 k - 1 comparisons. This is 
much more efficient than the "brute force" (i.e., complex) 
maximum likelihood decoding process which requires 
2 k x (2 m - 1) additions and 2 k - 1 comparisons. 
5 The Hadamard matrix H^ may be rewritten as a product 

of m matrices: 

H = E, x • x E 

ml m 

where E i is a (2 m x 2 m ) matrix with two nonzero (i.e., +1 
or -1) entries per column. In addition, for m'<m: 
10 F / = E, x • ■ ■ x E / 

ml m 

is a block diagonal matrix where the blocks are equal to 
H^. As discussed above, the fast Hadamard Transform (HT) 
is taken over the cosets for the 2 kU) versions of r' given 
by Equation (3) . An examination of the code words z 3 of 

15 the Reed-Muller subcode C {2) reveals that parts of the 
taken fast Hadamard Transform are needlessly recomputed 
during the recycling operation. Consequently, 
intermediate values of the Hadamard Transform are stored 
and used as needed to achieve significant savings in 

20 comput at ions . 

Reference is now made to FIGURE 2 wherein there is 
shown a flow diagram illustrating the decoding process 
implemented by the channel decoder 120 of FIGURE 1 for 
second order Reed-Muller codes. In step 200, the vector 

25 of the coded information is received. The process then 
precomputes and stores intermediate fast Hadamard 
Transform values in step 202. The decoder is then 
initialized in step 204 with j = -l. A determination is 
then made in step 206 as to whether all cosets have been 

30 exhausted (i.e., whether j = 2 k2 - 1) . If not, j is 
incremented in step 208. Then, in step 210, for the 
current value of j , an intermediate codeword guess is 
obtained over the coset C (1,j> of the subcode C ll) . The 
precomputed intermediate Fast Hadamard Transform values 

35 from step 202 are retrieved* and further processed to 
complete the Hadamard Transform to produce the 
intermediate codeword guess. The process then returns to 
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decision step 206. Once the process has cycled over all 
the cosets C a,i) of the subcode C aj , and found through 
completion of the Fast Hadamard Transform process an 
intermediate codeword guess at each cycle, the final 
5 codeword guess (comprising the best of the intermediate 
codeword guesses) is obtained therefrom in step 212. 
Decoding of the final guess is then performed in step 214. 

A more complete understanding of the present 
invention may be obtained now by examining the reduced 

10 complexity maximum likelihood decoding process for a 
specific second order Reed-Muller code. Take, for 
example, the second order (32, 16; 8) Reed-Muller code C 
(m=4, r=2) . One subcode C U) thereof comprises the first 
order (32, 6; 16) Reed-Muller code. Also, subcode C {2) 

15 comprises a (32, 10; 8) subcode of the second order Reed- 
Muller code. For the subcode C (2) , its code words are 
parsed into four blocks s, each of length eight. For 
8=1,2,3,4, let A s denote the set of distinct blocks of 
length eight taken by block s of code words in the subcode 

20 C (2) . Also, let a s denote the size of A g. Then, ap8, 
a 2 =64, a 3 =64, and a 4 =128. 

If a code word in the second order Reed-Muller code 
C is generated and transmitted, it is received as a 
receive vector r=(r 0 , . . . ,r n _ a ) . The vector r is then 

25 parsed into four blocks r l/ r 2 ,r 3 ,r 4 of length eight each. 

As an alternative, it is noted that parsing into eight 
blocks of length four each, or perhaps something in 
between, is also equally available. The same is true for 
two blocks of length sixteen each, or sixteen blocks of 
30 length two each. The choice of four blocks of length 
eight each yields the best compromise between computation 
and storage. For each r g/ a s Hadamard Transforms are 

computed and stored using the eight -by-eight matrix H 3 . 
Then, for each code word z^, the appropriate values of the 

35 oc s Hadamard Transforms are retrieved from memory, and used 
to complete the Hadamard Transform for multiplying the 
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retrieved values by E 4 x E 5 . In this regard, it noted 
that F 3 is a 32x32 block diagonal matrix where each block 
is equal to H 3 , and that H 5 =F 3 xE 4 x E 5 . Consequently, 
multiplying the intermediate values by E 4 x E 5 completes 
5 the Hadamard Transform of length thirty- two. 

Some simplifications are available. Looking first 
at the sets Ag, the additive inverse of each pattern in A 4 
is also in A 4 . Accordingly, the sixty- four included 
additive inverses may be deleted since their Hadamard 
10 Transforms need not be computed, giving A' s with a y 4 =64. 

Then, A' s = A s , and a' s = a s for s=l,2,3. It is next 
recognized that in some instances the pattern u in A' s 
coincides with a row of H 3 or its additive inverse (except 
for the first row) . Accordingly: 

(r g * «> H 3 = r s H 3 P (4) 

15 wherein: P is simply a permutation. 

Hence, there is no need to compute and store the Hadamard 
Transforms for such a pattern n in A' s/ giving A" s with 
a%=57 for s=2,3,4. Then, A" s = A' s , and a" ^ a' s for s=l. 
In summary, we then obtain new sets A" s with respective 

20 sizes of a- 2 =57, a" 3 =57, and a- 4 =57, and £ a s = 179. 

A full description of the operation of the channel 
decoder 120 may now be given. Let, u 3/i denote pattern 
number 1 in A" 3 . Then, let: 

For each code word z j of the Reed-Muller subcode C (2) , let 
25 Kj,s) denote the appropriate index in A" s pointing to the 
location of block s in z j; and let P(j,s) denote the 
permutation matrix to account, if needed, for Equation 
(4) . The maximum likelihood decoding of the received 
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vector by the channel decoder 120 for the second order 
(32 / 16;8) Reed-Muller code is then implemented by the 
following process: 

- compute t q1 , s=l,2,3,4, l-l #> .. r a", # according to 
5 Equation (5) , 

- 1=0, p=0 

- for j=0 : 2 k - 1 

- access stored values of ? |# s=l,2,3,4, 

- permute values using P(j,s), 8=1,2,3,4, if 
10 needed, 

- multiply by E 4 x E 5 to complete the 
Hadamard Transform of f , 

- find I / = arg max |r| , 

- if If^l > p, then 
15 - p = |* r ,|, 

- if r I > 0 then 1=1', otherwise I = I'+2\ 

- J-j, 

- end 

- end 

20 - guess y z + z a was the sent code word, 

- recover the information bits from I and J. 

Each H 3 requires 3x2 3 =24 additions, for a total of 
24x179=4,296 additions. Each E 4 x E 5 requires 2x32 
additions, and must be repeated 2 10 times, for a total of 

25 65,53 6 additions. Overall, this gives 69,832 additions. 
To store the intermediate Hadamard Transform values, 
179x8=1,432 memory units must be provided. Finally, there 
is a need for 2 16 -1 comparisons. This is much more 
efficient than the "brute force" (i.e., complex) maximum 

30 likelihood decoding process which requires 
2 16 x31=2, 031, 616 additions, or a ratio of twenty-nine 
times more additions than the reduced complexity 
implementation of the present invention. 

Utilizing the teachings herein, as well as the 

35 teachings of commonly assigned, co-pending U.S. 

Application for Patent Serial No. 08/768,530, filed 
December 18, 1996, by Ali S. Khayrallah, et al . , the 
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reduced complexity maximum likelihood decoding process of 
the present invention may be extended to shortened forms 
of second order Reed-Muller codes such as the second order 
(30 ,14; 8) shortened Reed-Muller code, or the second order 
5 (28 ,14; 6) shortened Reed-Muller code. Thus, to extend the 

efficient decoding process described above for use in 
connection with the decoding of shortened Reed-Muller 
codes, the generator matrix for the Reed-Muller code is 
modified to produce a modified generator matrix that is 

10 specific for and tailored to the decoding of each 
shortened code. The modified generator matrix is then 
efficiently implemented in the efficient decoding process 
to identify the best codeword for conversion to its 
corresponding information bits for output. In particular, 

15 the modified generator matrix comprises the Reed-Muller 
code generator matrix with specially chosen rows and 
columns deleted. 

The channel decoder 120 is preferably implemented as 
a specialized digital signal processor (DSP) or in an 

20 application specific integrated circuit (ASIC) . It will, 
of course, be understood that the channel decoder 120 may 
alternatively be implemented with the de-interleaver 118 
and source decoder 122 in combination or using discrete 
components and perhaps distributed processing. In either 

25 case, the de-interleaver 118, channel decoder 120 and 
source decoder 122 each perform and implement the 
functional operations previously described. 

Although preferred embodiments of the present 
invention have been illustrated in the accompanying 

30 Drawings and described in the foregoing Detailed 
Description, it will be understood that the invention is 
not limited to the embodiments disclosed, but is capable 
of numerous rearrangements, modifications and 
substitutions without departing from the spirit of the 

35 invention as set forth and defined by the following 
claims . 
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WHAT IS CLAIMED IS: 

1 . A method for maximum likelihood decoding of a 
received vector comprising second order Reed-Muller code 
encoded information bits, comprising the steps of: 
5 cycling over cosets of a subcode of a second order 

Reed-Muller code used to encode the information bits; 

performing a Hadamard Transformation on the received 
vector compared to each coset to obtain an intermediate 
codeword guess for each cycle; and 
10 recovering, following cycling through each of the 

cosets, the originally encoded information bits from a 
final codeword guess comprising a best one of the 
intermediate codeword guesses. 

15 2 . The method as in claim 1 wherein the Hadamard 

Transformation comprises a Fast Hadamard Transformation. 

3 . The method as in claim 1 further including the 
step of precomputing and saving a plurality of 

20 intermediate values for the Hadamard Transformation of the 
received vector. 

4 . The method as in claim 3 wherein the step of 
performing a Hadamard Transformation on the received 

25 vector further comprises the steps of: 

recalling necessary ones of the precomputed 
intermediate values at each cycle; and 

completing the Hadamard Transformation of the 
received vector using the recalled precomputed 
30 intermediate values to obtain the intermediate codeword 
guess with respect to each coset. 

5. The method as in claim 1 wherein the second 
order Reed-Muller code comprises a Kerdock code. 

35 
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6. The method as in claim 1 wherein the second 
order Reed-Muller code comprises a shortened form of the 
second order Reed-Muller code. 

5 7. The method as in claim 1 wherein the subcode is 

a first order Reed-Muller code. 

8 . A method for maximum likelihood decoding of a 
received vector comprising encoded information bits, 

10 comprising the steps of: 

cycling over cosets of a subcode of a code used to 
encode the information bits, the subcode comprising a 
first order Reed-Muller code; 

performing a Hadamard Transformation on the received 
15 vector compared to each coset to obtain an intermediate 
codeword guess for each cycle; and 

recovering, following cycling through each of the 
cosets, the originally encoded information bits from a 
final codeword guess comprising a best one of the 
20 intermediate codeword guesses. 

9. The method as in claim 8 wherein the Hadamard 
Transformation comprises a Fast Hadamard Transformation. 

25 10. The method as in claim 8 further including the 

step of precomputing and saving a plurality of 
intermediate values for the Hadamard Transformation of the 
received vector. 

30 11. The method as in claim 10 wherein the step of 

performing a Hadamard Transformation on the received 

vector further comprises the steps of: 

recalling necessary ones of the precomputed 

intermediate values at each cycle; and 
35 completing the Hadamard Transformation of the 

received vector using the recalled precomputed 
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intermediate values to obtain the intermediate codeword 
guess with respect to each coset . 

12. The method as in claim 8 wherein the code 
5 comprises a second order Reed-Muller code. 

13. The method as in claim 8 wherein the code 
comprises a Kerdock code. 



10 



14. The method as in claim 8 wherein the code 
comprises a shortened form of a second order Reed-Muller 
code . 
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